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Abstract 

Recent years have seen major advances in understanding the state of the intergalactic medium 
(IGM) at high redshift. Some aspects of this understanding are reviewed here. In particular, we 
discuss: (1) Different probes of IGM like Gunn-Peterson test, CMBR anisotropies, and neutral 
hydrogen emission from reionization, and (2) some models of reionization of the universe. 
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I. INDRODUCTION 



One of the outstanding issues in cosmology is to understand the reionization of the 
universe. Following recombination of primordial plasma at redshift z ~ 1000, the universe 



is mostly neutral with an ionization fraction ~ 10~ 4 (see e.g. [26], 2J, |2j|). The Jean's 
mass at recombination is ~ 1O 6 M . At z ^ 100, the plasma thermally decouples from 
CMBR and its temperature decreases adiabatically, T m ~ 1/a 2 , which leads to a further 
decrease in Jean's mass. During this 'dark and cold age' the density perturbations at scales 
above the Jean's scale can grow. Figure 1 shows the ionization and thermal history of the 
universe along with the evolution of Jeans' mass. This age comes to an end when the first 
structures can collapse and the light from these objects can reionize and reheat the universe 
(for a recent review see 2] and references therein) . Therefore the epoch of re-ionization 
holds important clues about the way first structures formed and can potentially distinguish 
between different models of structure formation. In the standard ACDM models the first 
structures to collapse would be just above the Jean's length (see e.g. j^J). Another crucial 
question in this regard is whether these structures could cool fast enough to form stars. 
Many of these issues will be discussed in this review. 

Different ongoing and potential probes of intergalactic medium can reveal the nature of 
the re-ionization of the universe. One of the most important and the oldest is the Gunn- 
Peterson test (see e.g. [3| and reference therein), which is very sensitive to the neutral 
fraction in the intergalactic medium. CMBR anisotropy measurements are another powerful 
and complementaryprobe, as they are sensitive to the ionized component of the intergalactic 
medium (see ^(|, Q] and references therein). In future, it might be possible to directly 
observe the first sources that re-ionize the universe. In addition, the transition from neutral 
;o ionised _ .night al80 be detected in nentra, hydrogen emission (see e., Q, Q, 

To sum up the observational status: Recent detection of temperature-polarization cross- 

n 

correlation in CMBR suggests that the redshift of reionization z re i on ~ 17 ± 4 [20]. Gunn- 
Peterson probes suggest that the universe is highly ionized upto z ~ 5, but might be 
making a transition from highly ionized to neutral for 5 £ z £ 6 0], ^j) . These two 

observations together throw open the interesting possibility that the universe went through 
two phases of re-ionization. 
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This article is divided into two parts. In the first part, probes that give a clue about the 
reionization epoch will be discussed. In the second part, we will discuss the nature of ionizing 
sources. Throughout this review we use the currently-favoured FRW model: spatially flat 



with Q m = 0.3 and Q\ = 0.7 



37j, |28|, with tt b h 2 = 0.02 aj, 0) and h = 0.7 



n n —i 



II. PROBES OF IONIZATION AT HIGH REDSHIFTS 
A. Gunn-Peterson effect 

High redshift sources should show absorption at frequencies close to Lyman-a (1216 A) 
line owing to scattering from IGM neutral hydrogen (HI). This test applied to high redshift 
quasars which generally have a strong Lyman-a emission line means that the blueward side 
of the Lyman-a line should show strong absorption as compared to the redward side (Gunn- 
Peterson {GP) test). The optical depth to the Lyman-a scattering from HI can be calculated 
(see e.g. I2J): 

-? ( h \ ( n b \ ( 0.3 \ 1/2 / 1 + z\ L5 

Here x m = n m /n u (n H ~ 0.92n b ) is the neutral fraction of hydrogen and z/ = v a j (1+z). Two 
points worth noting in Eq. (^) are: (1) Owing to the resonant nature of Lyman-a scattering, 
optical depth at observed frequency uq gives direct information about the neutral fraction 
at redshift (1 + z) = Va/fo-, and (2) more importantly, this test is extremely sensitive to the 
neutral fraction of hydrogen, even a neutral fraction as small as a part in hundred thousand 
can fully absorb light shortward of Lyman-a in the quasar spectrum. 

Since 1960s when the first high redshift quasars whose Lyman-a emission could be de- 
tected from ground-based telescopes were discovered, GP test has been applied to study the 
ionization of the universe at high redshifts. Till 2000, none of the observed quasars upto a 
redshift ~ 5 showed any GP absorption, which means the universe is ionized to better than 
a part in a million upto z — 5. The discovery of several quasars at redshifts above six by 
SDSS survey made it possible to apply GP tests at even higher redshifts. GP absorption has 
been detected in several quasars at redshifts z £ 5.7 (Q], j|, [l(|) . Observations suggest 
that the neutral fraction increases rapidly between z ~ 5.5 to z ~ 6.2 (see e.g. 10]). These 
observations might mean that the universe is becoming ionized owing to the formation of 
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FIG. 1: The temperatures of CMBR and matter are plotted alongside the evolution of Jeans mass. 
The solid curve is Jeans mass in units of 10 5 M . The dashed and dot-dashed curves correspond 
to matter and CMBR temperature (in Kelvin), respectively. 



the first structures at z ~ 6. However a straightforward interpretation of these results is not 
easy. Even for the best quasar spectrum and using the Lyman-/? line, which has a smaller 
oscillator strength, the GP optical depth r GP £ 25 (01)- Using Eq. (JTJ), this implies that 



10 4 i.e. the universe can be almost fully ionized. Use of semi-analytic models which 



take into account the dumpiness of the IGM give more stringent constraints on the neutral 
fraction and give x m ^ at z ~ 6 ^(|. Even though the observations are unable to con- 
clude that the universe made a transition from fully ionized to almost fully neutral between 
redshift of five and six, the ionized fraction certainly evolves very rapidly in this redshift 
range, and this can have important implications for the models of structure formation. 
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B. CMBR anisotropies 



CMBR anisotropies provide a complementary approach to the ionization history of the 
universe as compared to the GP test, as they are sensitive to the ionized component of the 
universe. The physics of CMBR temperature anisotropies at the last scattering surface and 
various data analysis issue are covered elsewhere in this volume (Subramanian, this volume). 
Here we shall discuss the implications of reionization on the CMBR temperature anisotropies. 
To highlight the effect of reionization on the polarization anisotropies we discuss more fully 
the CMBR polarization anisotropies. 

The universe recombines at z ~ 1000. Following recombination the ionized fraction (x e = 
n e/ n n) m the universe is ~ 10~ 4 for z ~ 100 (see e.g. The mean free path of the CMBR 

photons to Thompson scattering exceeds the local Hubble radius in the post-recombination 
era and therefore the universe is 'transparent' to the CMBR photons. Following reionization, 
the ionized fraction might reach nearly unity and a small fraction of CMBR photons might re- 
scatter again. An important quantity in studying CMBR quantity is visibility function which 
is the normalized probability that the photon scattered in a range z and z + dz. It is defined 
as: V(r] ,rj) = f exp(— r), here drj = dt/a is the conformal time and r = J® o x e n H aa t cdr]. 
In Figure 2 we show the visibility function for two models. For the model with early 
reionization, the visibility function get important contribution from redshift of reionization. 
One can define the optical depth to the reionization surface: r reion = J 2rcion x e n H a t cdt; this 
is the fraction of photons that re-scattered in the reionized universe for r re i on < 1. For 

^reion — 6, Trcion — 0.05 and T rc i on DC (l -f- Z re [ on ) ^ • 

The temperature and polarization anisotropies can be computed by solving the Boltz- 
mann equations for the photon distribution function. Equations appropriate for studying 
the effect of reionization for scalar perturbations, for a given wavenumber k and line of sight 
n, are (see e.g. [r|, Zalddariaga 1997, 4], ^]): 

A T + ik/jA T ~ f(fiv — At) 

A P + tkfiAp = f(n(/i)[A T2 + A P2 - A P0 ] - A P ) (2) 

Here /i = k.n, = dAx/drj, At = AT/T(k, fi, if). The polarization anisotropies 
Ap = AP/T(k, fi,r]) are in the plane perpendicular to the k vector, which is taken to 
be parallel to the z-axis, to exploit the axial symmetry of the problem. In this case, only one 
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FIG. 2: Visibility function, defined as dr / drj exp(— r)//^ 1 , is plotted for different models. The 
solid and the dashed curves are for the standard recombination and a model in which the universe 
reionizes at z = 17, respectively. 

Stokes parameter Q is non-zero and AP = QT. Also f = n b x e aTa, n(/i) = 0.5(1 — P 2 (A i ))- 
v is the electron velocity; we also assume curl-free velocity fields which allows us to ex- 
press: v.n = v.k. The angular moments of temperature and polarization anisotropies are: 
Ai = | j\ d^Pi{^)A{k 1 /i). Even though the temperature equation is only appropriate for 
studying the effect of reionization on anisotropies generated at last scattering surface, the 
polarization equation is exact and can also be used to study the generation of perturbations 
at the last scattering surface. 

To understand the essential physics of CMBR anisotropies in reionized models, we only 
solve for anisotropy evolution for one wavenumber k. The quantity measured by experiments 
is the two point function of this quantity summed over all the wave-numbers (for details see 
articles by Subramanian, this volume). 

We begin by studying the effect of reionization on the temperature anisotropies. Eq (J2J) 
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FIG. 3: The effect of reionization CMBR temperature anisotropies is shown. The main effect of 
reionization on temperature anisotropies is to generate secondary signal for I 20 and suppress 
primary anisotropies by a factor exp(— 2r re ion) for larger t. 

can be solved to give: 

A T (?7) = A T (r] rec ) exp(ikn(r] rec - 77)) exp(-r(r/ rcc , rj)) + n J drfv(rf)V(ri, rf) exp(ik/i(rf - r])) 

(3) 

The first term in the solution means the anisotropies generated at the last scattering surface 
are exponentially damped by reionization (the solution is only correct for scales smaller 
than the size of local horizon, i.e. k £ r]" 1 (z Teion ) . It is not reflected in the solution owing 
to dropping a term oc A^o in the temperature equation (for details see jlE^))- For example, 
if the universe reionized at z ~ 50 which gives r re i n — 1 -i-e. all the photons from the last 
scattering surface are re-scattered following reionizatio, this means that all anisotropies at 
scales smaller than the angular scale correspond to I ~ fjo/ rj(z IC i on ) — 10 are wiped out. As 
we shall see the second term in the solution doesn't contribute much to the anisotropies at 
small scales either. This means that for a reionization redshift ~ 50 no anisotropies should 
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be observed for i £ 10, which is in direct contradiction with observations (e.g. WMAP 
observations detect anisotropics for i ~ 600). Therefore the redshift of reionization should 
be small enough such that only a small fraction of CMBR photons are re-scattered. To 
compute the second term in Eq. ©, we can assume the visibility to be a normalized Gaussian 
with a width A^ reion , this gives 



This shows that CMBR anisotropics generated during the epoch of reionization owing to 
Doppler scattering off electrons are suppressed for k ^ l/A?7 re i on . The generated signal as 
expected is oc r re i n For realistic ionization histories A7/ reion — ^reion- We show in Figure 3, the 
effect of reionization on the temperature anisotropies. The net effect of reionization on the 
CMBR temperature anisotropies can be summarized as: (a) anisotropies at small scales are 
suppressed exponentially as exp(— 2r re i on ), (b) anisotropies at scales corresponding to I £ 10 
escape this exponential damping and also new anisotropies are generated at these scales. 
As seen in Figure 3, reionization causes a relative decrease in the small scale anisotropies. 
The minimum error in detecting the angular power spectrum of CMBR anisotropies at any 



£ is ACi ~ y 2/(2£ + l)Cg (Cosmic variance) and for all CMBR anisotropy experiments the 
error on Ce is dominated by the cosmic variance for £ ^ 300. The reionization signal is very 
difficult to detect owing to uncertainty in the overall normalization of the CMBR anisotropies 
which is compounded by cosmic variance. From WMAP data all the information about the 
reionization comes from the polarization signal. 

In second order in perturbation theory, reionization causes potentially detectable 
anisotropies for I ^ 1000; for i ^ 2000-3000 this signal can dominate the signal from primary 
anisotropies. The most important second order contribution to temperature anisotropies 

nnnn 

comes from Vishniac effect (for details see |4l(, also see [16j, [8j, |lja] for host of other second 
order effects). 

From Eq. @, it can be seen that the source of polarization anisotropy is the quadrupole 
of the temperature anisotropy. The generation of temperature and polarization anisotropies 
at the last scattering surface is generally an involved process and Eqs.(j2]) have to be solved 
numerically (see e.g. [s^)- However, several relevant assumptions allows one to get an 
analytic insight into the problem (Q|> Fhst simplification occurs because the photon- 

baryon plasma can be treated as tightly-coupled for most of the period during which the 



A T oc /!T reion f exp - ( kfiAr] reion /2 ) 2 



(4) 
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FIG. 4: Primary CMBR temperature (TT) and polarization (EE) anisotropies and their cross- 
correlation (TE) are shown; the absolute value of the temperature-polarization cross-correlation is 
plotted. 



recombination lasts. The tight coupling approximation is valid for scales corresponding to 
k £ Ij , where If is the (comoving) mean free path of the photons for Thompson scattering; 
If — 1 Mpc at the last scattering surface for a fully ionized universe, which means that tight 
coupling approximation is valid for most scales of interest. (The scale of interest for tight- 
coupling in not If but the scale of photon diffusion at recombination which for a fully ionized 
plasma is roughly 10 times larger than If] free-streaming of photons at recombination further 
increase this scale; for detailed discussions and implications of this for CMBR anisotropies 
see jlij). A second simplification, closely related to the first but physically distinct, occurs 
because the width of visibility function at the last scattering surface corresponds to scales 
£ 10 Mpc, and hence for studying physical processes at much larger scales the recombination 
can be treated as instantaneous (for caveats see 17|). Out interest here is in scales that are 
super- horizon at the last scattering surface (fc -1 ^ 100 Mpc, comoving), corresponding to 
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angular scales £ ^ 200. Therefore we will use tightly coupled approximation and not discuss 
the effects of photon diffusion. Using these approximations, adiabatic initial conditions and 
also the fact that the ratio of baryon to photon energy density p&/p 7 — 25Qi,h 2 <C 1 at the 
epoch of decoupling for acceptable models of primordial nucleosynthesis, the temperature 
and polarization anisotropies generated at the last scattering surface are: (for details see 

A T (k, n, rj) = Vrcc) cos(kc s r] rec ) exp[ikfi(r] rec - 77)] 

A P (k,n,T)) = 0.17(1 - yU 2 )A;c s A?7 rcc $(k, r/ rec ) sm(kc s i] rcc ) exp[zfc/i(r/ rcc - rj)] (5) 

Here $(k, ?7 rec ) is the Fourier component of the Newtonian potential at the last scattering 
surface; A?7 rec ~ lOMpc is the comoving width of LSS; c s ~ l/v3 is the sound velocity in 
the coupled photon-baryon fluid. These large scale solutions show that: (a) temperature and 
polarization anisotropies are correlated with each other, (b) the amplitude of polarization 
anisotropies is suppressed by a factor kc s Ar] Tec ; e.g. for k ~ 0.01 Mpc -1 this factor is roughly 
1/10. This suppression is owing to the fact that polarization anisotropies oc A^ 2 - In the 
strict tight coupling approximation, only the temperature monopole and dipole are non-zero 
in the comoving frame. A small quadrupole is generated owing to free streaming of photons 
before they scatter for the last time; and this quadrupole is suppressed with respect to the 
monopole and dipole, which constitute the primary sources of temperature anisotropies, by 
the factor ~ kc s Ar] Tec . As we shall see below, in the reionized models this suppression is 
absent. In Figure 4 we show the CMBR temperature and polarization anisotropies, generated 
at the last scattering surface. 

In the reionized models, a fraction CMBR photons scatter again after the epoch of reion- 
ization. As argued above, this fraction is small as many of the features of the primary 
anisotropies have already been detected. We have already discussed the effect of reioniza- 
tion on temperature anisotropies and argued that the reionization signal is very difficult 
to discern from temperature anisotropies alone. Like temperature anisotropies, one of the 
effects of reionization would be to wipe out polarization anisotropies generated at the last 
scattering surface for scales £ ?7reion- From Eq. (j2J), the generation of new polarization 
anisotropies will be proportional to the value of At2 at the epoch of reionization. Eq. (J2J) 
can simplified further by dropping terms of polarization monopole and quadrupole of po- 
larization anisotropies in the RHS of the equation, as these terms are negligible compared 
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FIG. 5: The effect of reionization on the temperature-polarization cross-correlation is shown. Also 
shown are the WMAP data points with one-sigma error bars [2o| 

to the temperature quadrupole (see below). Eq. (J2J) then give the following solution for the 
generation of polarization anisotropies following reionization: 



dr]'Yl(fi)A T2 (ri')V(r], rj') exp(ikp,(r}' - 77) ) 



(6) 



Most of the contribution to the integral will come from close to the reionization epoch 
(typically ~ A?7 rc i on ~ ?7reion, the width of visibility function if it is approximated as a 
Gaussian as above). The amplitude of the contribution is proportional to the temperature 
quadrupole at the epoch of reionization. As argued above the temperature quadrupole 
is suppressed at the last scattering surface owing to tight coupling approximation that 
holds at that epoch. Following recombination the photons free stream which allows to 
temperature quadrupole to increase. From Eq. (JHJ, we can get the temperature quadrupole 
at the reionization epoch from taking the angular moment of the equation (for details see 
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Ay (r] vcion ) = ^roc) cos(kc s 7] rcc )j 2 [k(r] ieion - r} Iec )) (7) 

Using the fact that the maximum value of j2{x) ~ 0.3, and comparing with the polarization 
anisotropies generated at the last scattering surface, Eq. (JHJ), we can see that (''/reion) 
doesn't suffer the k dependent suppression, and is of the order of temperature anisotropies 
at the last scattering surface. At k -C 0.1, it can be several orders of magnitude more 
than the polarization anisotropies generated at the last scattering surface. The polarization 
anisotropies at the reionization epoch peak at the characteristic scale k ~ 2/i] rcion , which 
corresponds to an angular scale I ~ kr] ~ 5-10 for z re ion — 15-50. 

In Figure 5 we show the effect of reionization on the temperature-polarization cross- 
correlation power spectrum; also shown are the observations of WMAP. Figure 5 shows 
that the enhancement of power at £ < 10 cannot be explained within the framework of 
no reionization models. Also shown are the predictions of a model in which the universe 
reionized at z ~ 5.5, which the GP observations discussed above might be suggesting and 
the best fit model to the WMAP observations. The best fit model requires r re ion — 0.15, 
which implies the epoch of reionization corresponds to z re ion — 15. 



C. High Redshift HI 

Another possible probe of the reionization epoch is to observe the neutral component of 
hydrogen thorough the epoch of reionization. The neutral fraction of hydrogen changes from 
near unity to zero during the epoch of reionization. This change can potentially be observed 
using the hyperfine transition of the hydrogen atom at v+ = 1420 MHz. The quantity of 
interest here is the spin temperature of hydrogen defined as: 

^ = 3exp(-T,/T s ) (8) 

m 

Here n 2 and ni are the populations of the hyperfine states. T± = hv^jk = 0.06 K. As the 
only radio source at the high redshift is CMBR, HI in hyperfine transition can seen against 
the CMBR in emission or absorption. The observed quantity then is the deviation of CMBR 
from a black body at radio frequencies. The observed difference is: 

ATcmbr = — t H i(Tcmbr — T 8 ) (9) 
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FIG. 6: ATcmbr (Eq. ©) is shown as a function of the observed frequency for an ionization 
history in which the universe becomes fully ionized at z = 15 

Here r m = a u NmT^/T s , with the HI column density iV H i = / n m d£; o v = c 2 A 2 i0jy/(47rz/ Jr ); 
A21 — 1.8 x 10 _15 sec _1 and <p u is the line response function. The spin temperature is 



determined from detailed balancing between various processes that can alter the relative 
populations of the two levels jll|, [l^ : 

^cmbr + VcXk + Va^a 



s ~ l + y c + y« 1 ' 

Here y c oc ri H and y a oc n a , with n a being the number density of Lyman-a photons cor- 
respond, respectively, to relative probabilities with which collisions between atoms and the 
presence of Lyman-a photons determine the level populations; Tk is the matter temperature. 
In the pre-reionization era, there are no Lyman-a photons, and therefore y a = 0. 

From Eq. (jHJ) it is clear that HI can be observed in either absorption or emission against 
CMBR depending on whether T s is less than or exceeds Tcmbr- At z ~ 1000, Tcmbr = Tk 
and it follows from Eq. (fTUj) that AT C mbr = 0. As seen in Figure 1, for z £ 100, T K < T C mbr 
and therefore HI can be observed in absorption against CMBR if y c £ 1. During reionization, 
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y a term can become important and owing to thermalization at frequencies close to Lyman- 
a, T a ~ Tk [l^- The temperature of the medium can also exceed Tcmbr from X-ray 



and Lyman-a heating (for details see |2j|, |33j) which implies that AT CM br > 0. AT CM br 



approaches zero as the reionization is completed. In Figure 6, we show ATcmbr as a function 
of the observed frequency for an ionization history in which the universe becomes fully ionized 



at z = 15 (more details in next section and 



3)- 



If the universe re-ionized at z ~ 15, then 



ATcmbr — 0.05 K at frequencies v ~ 50-80 MHz. In addition there is a signal from pre- 
reionization epoch with ATcmbr — — 0.05K at v ~ 30 MHz. In addition to the average 
signal there will be fluctuations in the temperature difference owing to fluctuations in HI 
density from primordial density perturbations and also from the patchiness of reionization 
(for details see e.g. [s^). Currently the prospects of detecting this signal are being studied 
by using both single dish and interferometric experiments at low radio frequencies ([35], see 
also [23). 

III. RE-IONIZATION OF THE UNIVERSE AT HIGH REDSHIFT 

In the previous section we discussed various probes of the high redshift universe and the 
epoch of reionization. In this section we take up the issue whether it is possible to explain 
the observed ionization structure of the universe within the framework of currently-favoured 
ACDM models of formation of structures. In these models, the observed structures in the 
universe grew from gravitational instability of density fluctuations which originated dur- 



ing inflationary epoch in the very early universe (see .e.g 26|, |24l. l25| ) . The gravitational 
collapse of these structures might either set off star-formation (alternatively some material 
might end up in black holes which by accreting more matter will radiate with harder spec- 
trum than first star-forming galaxies, see e.g. j29() which will emit UV light and ionize the 
IGM. The process of re-ionization of the universe is generally quite complicated and not 
well understood. However it is possible to study it within the framework of simple models 
which might give important clues about the details of this process. Important ingredients 
of this problem are: (a) Halo population at high redshift, (b) molecular and atomic cooling 
in Haloes, (c) Initial Mass Function of stars and star formation rate, (d) Escape fraction of 
UV photons from Haloes, (e) dumpiness of the IGM. Of these the most uncertain are (c) 
and (d) and have to be modelled using simple parameterized model. 
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Halo Population: The number density of dark matter haloes per unit mass M at any 



redshift can be obtained from the Press-Schechter method (see .e.g 
by: 



124 



251] ). It is given 
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dM 
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-5 c {zf/{2a\M)) 



(11) 



a 2 (M) 

Here a(M) is the mass dispersion filtered at the scale corresponding to mass M in the linear 
theory; a(M) for length scale corresponding to 8/i~ x Mpc is ~ 0.9 [37]; 5 c (z) ~ 1.7D(0)/D(z), 
with D(z) being the solution of growing mode in linear theory; in sCDM model D(0) / 'D(z) = 
(1 + z) (for details see 24, (25], j^]). The first baryonic structures to collapse would have 
masses exceeding the Jeans mass of the IGM, whose evolution is shown in Figure 1. From 
Figure 1 it is seen that the first structures might have masses ~ 10 4 M (for non-linear 
extension of the concept of Jeans mass see e.g. [2J; the modified mass has values similar 
to the one obtained using linear theory). The collapsed fraction of all structures can be 
calculated from Eq. (JTTj) and is ~ 2 x 10~ 3 at z simeq20. However to form stars, baryons 
need to cool sufficiently rapidly in the dark matter haloes to collapse to density higher than 
initial virialized density. 

Atomic and molecular cooling: Assuming a spherical, top-hat collapse (see e.g. [24! I25I]) 
the density of the collapsed structure is ~ 170 times the background density at that redshift. 
If a virial equilibrium is reached inside the halo, the baryon temperature is raised to the 
virial temperature given by (e.g. j3, 13) : 
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Here molecular weight ji = 1.22 for a fully neutral halo (haloes with masses £ 1O 8 M at 
z — 20); /x = 0.57 for a fully ionized halo. An important criterion is whether baryons can 
cool rapidly enough so that they collapse to higher densities, fragment and form stars. In 
the primordial gas, the cooling in haloes with virial temperature £ 10 4 is dominated by 
atomic hydrogen and singly- ionized helium. But the smaller haloes can only cool further by 
molecular hydrogen. A small fraction ~ 10~ 6 of molecular hydrogen is formed at z ~ 1000 
in the IGM following recombination of the universe (see e.g. |26j]). Such a small fraction 
doesn't suffice to cause rapid enough cooling. However the collapse of halo can cause the 
formation of molecular hydrogen, resulting in a molecular fraction ~ 5 x 10~ 4 by the time 
the halo virializes (see e.g. 39j and reference therein). In Figure 7, we show the cooling 
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FIG. 7: Cooling time is shown against the mass of the virialized halos at z = 20. The solid line 
is from atomic cooling of the primordial gas containing only hydrogen and helium. The cooling 
processes include: line cooling from hydrogen and helium, free-free emission, and recombination. 
The dashed curve is for cooling from molecular hydrogen for a molecular fraction 5 x 10 -4 . The 
dotted line is the inverse of local expansion rate H~ 1 {z) and the dot-dashed line is the dynamical 



time Gp n 



time (defined as t coo i = kTx/E, where E is the rate at which the halo loses energy from 
various processes) for haloes of different masses at z — 20. The criterion for runaway cooling 
resulting in fragmentation and star formation is that cooling time be less than the dynamical 
time (td = 1/y/Gp) of the halo. From Figure 7 it is seen than haloes of masses ~ 10 6 M 
can collapse to form stars at z ~ 20. This result is only approximate as the haloes will not 
have constant density. For instance, the more dense central parts of the halo can form both 

n 

more molecular hydrogen and cool more rapidly. Bromm et al. p] simulated the collapse of 
haloes of masses ~ 10 6 M Q with substructure. They concluded that the halo fragments into 
many clumps with typical masses 1O 2 -1O 3 M ; their analysis suggests that this mass scale 
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FIG. 8: 

dashed 
iV 7 (0) = 



20 
Redshift 

Ionization history is shown for several values of iV 7 (0) and C (see text for details) 
and dot-dashed curves are: {iV 7 (0) = lO^sec" 1 , C = 1; 7V 7 (0) = 2 x lO^sec" 1 , 
= 5 x 10 50 sec _1 , C = 5}, respectively. 
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. Solid, 
C = 2; 



is a result of molecular hydrogen chemistry and therefore should correspond to the masses 
of the first stars. Abel et al. |l| reached similar conclusions from their simulations. 

From information about the halo population and cooling arguments it is possible to 
speculate on the ionization history of the universe from photo-ionization. The main un- 
certainty is the hydrogen-ionizing luminosity (and its evolution) of each halo, which has to 
be parameterized. Assuming that halo of mass M emits isotropically the hydrogen-ionizing 
luminosity iV 7 (in photons sec -1 ), the radius of ionizing sphere around the source will satisfy 
the equation (Stromgren Sphere, see e.g. js^J, j^ll): 



d 4-HR 
dt 



(A> 7 - An/ZFPctBCnlxt 



(13) 



(N 7 + A-nR 2 x m n b ) 

Here C = (nf)/ (rib) 2 is the clumping factor of the IGM. Using Eq. ([lip, the fraction of the 
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47r c z r (In 
fUz) = -j o dz' j dM—{M,z>)R\M,z,z') (14) 



universe that is ionized as a given redshift is (see e.g. 

4-7T f z . , r dn 



Further assuming that the luminosity of the source is oc M, the ionized fraction can be 
calculated in terms of the evolution of the photon luminosity of a single halo of some fiducial 
mass and the evolution of the clumping factor. For simplicity we assume the clumping 
factor to have a constant value between one and five, and take the luminosity evolution of 
a halo to have the form of a typical star-burst galaxy: iV 7 (t) = iV 7 (0) exp(— 1/10 7 years). In 
Figure 8 we show several ionization histories for different values of A^,(0) for M = 5 x 10 7 M 
and clumping factor C. In Figure 8 we only plot ionization fraction upto z ~ 9 at which 
f ion ~ 1; further evolution keeps the universe fully ionized upto the present. These ionization 
histories are consistent with WMAP observations. (If the universe becomes fully ionized at 
z ~ 8 if will remain ionized upto the present even in the absence of ionizing sources as the 
typical recombination time 1/ (asCrib) already exceeds the age of the universe at this epoch 
for C £j 2.) It appears that iV 7 (0) ~ 10 50 is required to ionize the universe early enough 
to satisfy WMAP observations. This is just three orders of magnitude below the photon 
luminosity of a typical star-burst galaxy (see .e.g.j^lj]). Alternatively one can infer that 
the efficiency of the first star formation was very high (see e.g. Q], p). There are other 
uncertainties like feedback from supernova, and photo-dissociation of molecular hydrogen 
which indicate this estimate is a lower limit (see e.g. Q|, [l^, |?|). Alternatively it is 
possible that the collapsed fraction of the universe far exceeded the value given by ACDM 
models and was caused by some other physical process like tangled magnetic fields 3^. 

Is it possible to simultaneously satisfy the WMAP and the GP observations? It appears 
difficult to achieve it unless the efficiency of star formation decreased rapidly from z ~ 20 
to z ~ 6 (for other possible scenarios and details see e.g. Q], and references therein). 



IV. DISCUSSION 



The ionization history of the universe at high redshift as inferred by WMAP and GP 
observations is quite complicated: the universe reionized at z ~ 12-15; the neutral fraction 
increased to better than one part in a thousand for 5.5 < z < 6.5; and for z £j 5, the neutral 
fraction dropped to one part in a million. Theoretical analyses based on models of structure 
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formation are not sophisticated enough to understand this ionization history. Future obser- 
vations however might throw light on these observations. Firstly, future telescopes might 
be able to detect the sources of reionization upto z ~ 20 Future CMBR experiment 

°lanck also is sensitive enough to disentangle the effects of complicated ionization histories 
3- 
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